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Abstract 

Using a Hubbard-Stratonovich transformation coupled with Fourier path 
integral methods, expressions are derived for the numerical evaluation of the 
microcanonical density of states for quantum particles obeying Boltzmann 
statistics. A numerical algorithm is suggested to evaluate the quantum density 
of states and illustrated on a one-dimensional model system. 



Typeset using REVTgX 



Over the past decade there has been considerable progress in the development of path 
integral [|I|] approaches to computational quantum statistical mechanics [0,0. With few 
exceptions |3J] the application of path integral methods have been restricted to simulations 
in the canonical ensemble. Canonical simulations are amenable to path integral treatments 
because the canonical density matrix is formally identical to the quantum propagator in 
imaginary time. The purpose of this note is to show that path integral methods can provide 
an algorithmic basis for microcanonical simulations as well. 

We let Q(E)dE represent the number of energy states between E and E + dE. Through- 
out this paper we will suppress the dependence of the density of states on the number of 
particles in the system N and the system volume V for notational convenience. For N indis- 
tinguishable particles obeying Boltzmann statistics, we begin with the Fourier path integral 
expression for the canonical partition function [|J 
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In Eq. ([!]) we have suppressed the dependence of the partition function on N and V, 
(3 = 1/ksT with T the temperature and ks the Boltzmann constant, r is a collective 3N- 
dimensional vector representing the coordinates of all the particles in the system, 

m is the mass of the constituent particles, a^j is Fourier coefficient k for particle j, and for 
the potential energy $(r) we define 

< $ > a = f 1 du$[r(u)} (3) 

as the average of the potential energy over a particular path. As discussed elsewhere |J), 
each path is parameterized by a Fourier sine series. In practice a finite number, k max , of 
Fourier coefficients are actually included in Eq.(|l]). 

The connection between the canonical partition function and the microcanonical density 
of states is given by 
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Q((3) = J dEn(E)e-? E (4) 

which is just a Laplace transform relation. Equation (|j) can be inverted using the Mellin 
inversion integral 

Q(E) = 7t~t [ d^ E Q{(3) (5) 
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where the subscript C on the integral in Eq. (|5D implies integration along a contour to the 
right and parallel to the imaginary axis in the complex plane. The inversion in the present 
application is not simple, because of the inverse /3-dependence in the exponent of Eq. ([[]). 
The /3-dependence in Eq.(|l]) can be altered using a Hubbard- Stratonovich transformation ||. 
Writing 

«p(- E = II (^f) ' I d <t> ex P(" E <A, + »V2 E »«^) ( 6 ) 
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the expression for the partition function becomes 
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where 



°h = P°lr (8) 

At the expense of the introduction of the auxiliary variables 4> , the /5-dependence of the 
exponent in the expression for the canonical partition function appears entirely in the numer- 
ator. Using Eq. (|7|) the canonical partition function can be introduced in the Mellin inversion 
integral, and the integration with respect to (3 can then be evaluated analytically || 
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where T(x) is the gamma function and Q(x) is the step function. Equation is the principal 
result of this note. It is easy to verify that Eq.(|) reduces to the proper free particle and 
classical limits. 

A possible approach to evaluate Eq.@ numerically is based on the existence of 
methods to compute the classical density of states. We then correct the classical density of 
states by calculating the ratio of the quantum to classical density. If we specialize Eq. (f|) to 
a one-dimensional system and examine the classical limit, it is easy to show that the ratio 
of the quantum to classical density of states can be written 

n d (E) \y/2n) <5( a )5(<f>)> [ ' 

where the averages are taken with respect to the weight function 

In Figure 1 we present the quantum density of states for a Morse oscillator = D e {l — 
exp[— a(x — x e )]} 2 . The calculations in Figure 1 were evaluated from Eq. ([370]) and consisted 
of one hundred million Monte Carlo points with k max = 2. The classical density of states 
was determined from the classical limit of Eq. (|9]) and was evaluated using Mathematica || . 
For the potential parameters we have chosen a = 1, x e = 1, m = 1822.83 au, and we have 
chosen D e so that the oscillator frequency cu = (2a 2 D e /m) 1,/2 equals .006 au. The delta 
functions in Eq.flToD were represented by gaussians with standard deviations of 0.5 au. The 
path average of the potential energy for the Morse oscillator was evaluated with trapezoid 
rule quadrature with 16 quadrature points. The peaks in the density of states seen in Figure 
1 are in good agreement with the analytic eigenvalues at E/hu = .488, 1.40 and 2.21. The 
finite widths of the peaks are a result of the gaussian representation of the delta functions 
and the finite number of Fourier coefficients used. These preliminary results are encouraging, 
and the investigation of several methods to evaluate Eq.(|9|) for more complex systems is in 
progress. 
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FIGURE CAPTIONS 



1. The quantum density of states in atomic units for a one-dimensional Morse oscilla- 
tor with parameters defined in the text. The error bars are at the double standard 
deviation level. The maxima in the density of states are in good agreement with the 
analytic eigenvalues at E/Tiw = .488,1.40 and 2.21. 
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